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ABSTRACT 

Several  synoptic  problems,  previously  solved  by  graphical  overlays,  can 
be  solved  with  a  programmable  calculator.  The  HP-41CV  hand-held  calculator  is 
used  to  determine  surface  winds,  the  probability  of  average  rainfall  for  a 
regional  area  (Monterey,  California),  contrail  formation,  and  the  lifting 
condensation  level  from  surface  data. 


IV 


I.   Introduction 

This  report  describes  four  programs  that  use  the  HP-41CV  calculator's 
programming  capability  in  dealing  with  interesting  meteorological  problems, 
including  determining  surface  winds,  the  probability  of  average  rainfall  for  a 
regional  area  (Monterey,  California),  contrail  formation  and  the  lifting 
condensation  level  from  surface  data. 

The  programs  were  prepared  to  use  the  hand-held  calculator  in  solving 
synoptic  problems.  The  first  problem  is  the  estimation  of  surface  wind  at  a 
point  from  a  regional  sea-level  pressure  analysis  or  forecast.  This  problem 
arose  in  the  synoptic  laboratory  where  wind  forecasts  were  prepared  with 
numerical  weather  prediction  guidance  (facsimile  sea-level  pressure  maps). 
The  immediate  solution  is  a  geostrophic  wind  scale  applied  to  the  facsimile 
chart.  The  HP  program  provides  an  alternative  solution.  The  program 
calculates  the  geostrophic  and  gradient  wind,  as  well  as  the  estimated  surface 
wind,  using  Ekman  boundary  layer  assumptions,  by  simply  entering  latitudes, 
longitudes  and  pressure  readings  from  a  sea-level  pressure  chart. 

The  second  problem  involves  preparing  a  statistical  precipitation 
forecast  for  a  set  of  graphical  tables.  Using  a  tutorial  program  in  a  hand- 
held calculator,  the  need  for  the  graphical  tables  is  eliminated.  The  program 
determines  the  probability  of  average  rainfall  over  the  Monterey,  California 
area  for  the  period  0800  PST  forecast  day  to  0800  PST  of  the  following  day, 
for  the  wet-season  months,  November  through  April.  This  procedure  is  based  on 
the  Renard  (1972)  graphical  objective  technique  for  forecasting  24-h  rainfall 
at  Monterey,  California.  From  the  data  entries  of  500  mb  heights  used  to 
calculate  geostrophic  relative  vorticity  at  Monterey  and  a  point  8°  latitude 
upstream  from  Monterey,  as  well  as  sea-level  pressure  data,  the  program  uses 
a  series  of  pre-established  data  tables  to  arrive  at  a  probability  of  24-h 
rainfall  for  Monterey,  California. 


The  third  problem  involves  forecasting  contrails  from  upper-level 
sounding  data.  It  was  desired  to  have  the  HP  make  the  decision  from  input 
sounding  data  as  an  alternative  to  plotting  the  entire  sounding  and  using  a 
graphical  overlay.  This  program  determines  contrail  formation  at  a  given 
pressure  level  based  on  the  curves  found  in  the  Contrails  Forecasting  Manual 
of  Chief  of  Naval  Operations  (1964).  A  least  squares  approximation  was  made  to 
the  curves  involving  an  expression  for  critical  temperature  which  is  a 
function  of  temperature,  pressure  and  relative  humidity.  If  the  relative 
humidity  is  estimated  (dew-point  temperature  is  unknown),  there  is  a  series  of 
criteria  to  determine  whether  or  not  contrails  will  occur. 

The  final  problem  is  to  determine  the  lifting  condensation  level  (LCL) 
without  using  a  thermodynamic  diagram.  Inputing  values  of  surface 
temperature,  surface  pressure  and  specific  humidity  to  the  HP  results  in  a 
lifting  condensation  level.  The  LCL  is  the  height  at  which  a  parcel,  when 
lifted  dry-adiabatically,  becomes  saturated.  The  program  uses  the  Clausius- 
Clapeyron  equation  to  arrive  at  a  final  expression  for  the  LCL  that  is 
determined  by  values  of  surface  temperature,  surface  pressure  and  specific 
humidity. 


II.   Computation  of  Surface  Winds 

A.  Objectives 

This  program  will  give  the  forecaster,  either  ship-board  or  at  a  regional 
center,  an  easy  and  accurate  method  of  computing  the  geostrophic,  gradient 
and  estimated  surface  wind  from  a  MSL  pressure  analysis  or  model  .forecast. 
This  method  is  more  complete  than  a  graphical  geostrophic  wind  scale  overlay. 

B.  Principles 

The  geostrophic  wind  speed  and  direction  are  obtained  from  the  east-west 
(ug)  and  north-south  (vg)  components. 


ug  =  - 


_1   3P 

fp  3y 


|vg|=  Ju< 


g     +  v2g 


1   3p 

V*  '  fp  3x 


Vg 

t  an  6  g  =  — 

Ug 


where  ug  and  vg  are  geostrophic  wind  components, 

f  =  2  Q   sin  <t>  ;   is  latitude  of  forecast  point, 
p  =  1.3  kg/m3 ;  typical  sea  level  density  value, 
9x,3y  =  1°  lat  =  111  km 


Figure  1 

Geostrophic  wind  fVg )  computed  by 
the  finite  differencing  scheme 
with  1°  latitude  grid  spacing. 


The  user  enters  the  latitude  and  longitude  of  the  point  (Northern 
Hemisphere)  where  a  wind  is  desired.  Next,  the  user  enters  four  pressures 
that  are  1°  latitude  north,  south,  east  and  west  of  the  center  point  used  in 
the  finite  differencing  technique: 

dp   .   P@-  P(?)   3P    pQ-  P®      .„ 

3y  "   222ta    :  8x  =  "mta    '    (See  FlgUre  1*) 

The  gradient  wind  speed  is  obtained  from  the  geostrophic  wind  speed  and 
the  radius  of  curvature  (R) . 

■   .  [VgJ where  R  >  0    cyclonic 

'  *r'     .5  +  ^25  +  Vg/fR  R  <  0    anticyclonic 

The  user  enters  a  radius  of  curvature  (in  kilometers)  which  is  representative 

of  the  isobaric  curvature  in  the  forecast  area.   (See  Figure  1.) 

In  a  cyclonic  situation  (R  >  0),  the  magnitude  of  the  gradient  wind  is 
less  than  the  magnitude  of  the  geostrophic  wind  (subgeostrophic) ;  if  R  <  0 
(anticyclonic),  the  gradient  wind  speed  is  greater  than  the  geostrophic 
(supergeostrophic) .  If  the  curvature  is  zero  (R  =  <» )  the  magnitude  of  the 
gradient  and  geostrophic  winds  is  the  same.  Regardless  of  the  curvature,  the 
gradient  wind  is  from  the  same  direction  as  the  geostrophic  wind. 

The  estimated  surface  wind  speed  and  direction  are  obtained  from  the 
gradient  speed  and  direction  using  Ekman  layer  assumptions.  Over  the  land, 
the  gradient  wind  speed  is  multiplied  by  a  factor  of  0.6  which  is  an 
approximate  amount  that  the  speed  is  reduced  by  friction.  The  resulting  speed 
is  an  estimated  surface  wind  speed.  The  actual  direction  is  rotated  counter- 
clockwise 20°  from  the  original  direction.  Over  the  sea,  the  gradient  wind 
speed  is  multiplied  by  .81  and  direction  is  rotated  by  10°.  The  oceanic  and 
land  frictional  effects  follow  the  observational  summary  of  Sheppard  (1969) 
and  the  discussion  of  Petterssen  (1956)  and  Haltiner  -and  Martin  (1957). 
Sheppard  ('1969)  and  Mendenhall  (1967)  indicate  10°  is  a  reasonable  estimate  of 


frictional  veering  over  the  ocean  from  weather  ship  data.  Over  land,  the 
frictional  effects  are  naturally  heavily  dependent  upon  surface  roughness  so 
the  20°  veering  is  only  a  general  estimate. 

C.   User's  Guide 

1)   Make  sure  the  calculator  and  printer  (NOR  mode)  are  on. 

SIZE 


2) 


6) 


10 
11 
12 
13 
14 
15 
16 


XEQ   alpha 


alpha 


will  prompt  SIZE  

3)  Type  100. 

4)  Enter  magnetic  cards  as  labelled  1-9  -  - 

5)  After  the  program   (WIND)  is  "WORKING"  and  entered,  it  is  ready  to  be 
executed. 

WIND 


XEQ   alpha 


alpha 


7)  The  prompt  CTR  LAT?  will  appear. 

8)  Press  the   R/S  button  to  continue. 


(Do  this  each  time  you  answer  a  prompt.) 
9)   The  prompt  CTR  LON?  will  appear. 

Enter  the  longitude  of  the  center  forecast  point. 
Press 


R/S 


Answer  the  prompts  for  the  four  pressures  (mb). 
(Remember  to  do  this  each  time.) 


Press 


R/S 


The  geostrophic  wind  speed  (m/s)  and  direction  will  be  displayed. 

Enter  a  radius  of  curvature  (m)  after  prompt  RADIUS? 

The  gradient  wind  speed  will  be  displayed. 

The  prompt  LAND  or  SEA?  will  appear. 

If  the  forecast  region  is  over  land,  enter  "LAND";   otherwise,   enter 

"SEA". 


17)   The  estimated  surface  wind  speed  and 
the  program  is  done. 


direction  will  be  displayed  and 


18)   To  run  the  program  again:   XEQ   alpha   WIND   alpha 


D.   Examples 

To  illustrate  the  use  of  the  program  the  surface  wind  will  be  estimated 
for  points  A  and  B  in  Figure  2.  The  observed  surface  wind  speed  and  direction 
for  forecast  point  A  at  0000Z  on  26  November  1981  is  10  m/s  from  315° .  For 
point  B,  the  surface  wind  speed  is  20  m/s  from  295°. 

Execution  of  the  "WIND"  program  on  the  HP-41CV  calculator  results  in  a 
surface  wind  speed  of  14  m/s  from  a  northwest  direction  of  319°  for  forecast 
point  A.  Entering  the  latitude/longitude  position  (43N  130  W)  along  with  four 
pressures  1°  latitude  north,  south,  east  and  west  (1007  mb,  1012  mb,  1009  mb, 
1012mb  respectively)  results  in  a  geostrophic  wind  of  20  m/s  from  329° .  A 
radius  of  curvature  of  1,100  km  allows  a  gradient  wind  to  be  computed  (18 
m/s).  Finally,  since  forecast  point  A  is  an  ocean  location,  the  oceanic 
Ekman  layer  boundary  assumptions  are  applied.  The  resulting  estimated  surface 
wind  is  14  m/s  from  319° .  This  estimated  wind  is  in  agreement  with  the 
observed  real   wind  speed  (10  m/s)  and  direction  (315°). 

For  forecast  point  B  (46N  135W)  the  four  pressures  (1008  mb,  1011  mb, 
1006.5  mb  and  1013  mb)  result  in  a  geostrophic  wind  of  24  m/s  from  295° . 
From  the  radius  of  curvature  (1,300  km),  a  gradient  wind  of  21  m/s  is 
calculated.  The  final  estimated  wind  speed  for  this  ocean  location  is  17  m/s 
from  285°  which  compares  with  the  observed  real  wind  speed  (20  m/s)  and 
direction  (295°). 


Figure  2 

Surface  conditions   observed  at   0000Z  26  November,   1981   for 
illustration  of  surface  wind  program. 
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I.       Program 


Table  I 
Surface  Winds  Program. 

=  ♦LBL  "WIND*  /  li?!  PTmh  151  sTQ  18 

FIX  0  1*3  182  186  152  GTO  "GRPD" 

"CTf\  LPT?"  -  '  Iw3  +  153*LBL  "EPST- 

*4  KKUflPT  54  RCL  «3S  t@4  c,tq   js  j^,   .^,h 


is 


95  STO  8t  5j  /  iQ5  qtq  "G?aD*  155 

06  "CTR  LON?"  5*1  STO  18  106+L8L  "^SH"  156  GTO  "GRPD" 

ft?  PRQHFT  57  RCL  81  jfi?  m  157*lBL  -GRPD" 

m  STO  82  58  EHTERt  195  ENTER!  158  "GEO.  DIR  PRO»t 

99  -pS'S  HO.  LflT?"  59  Kt2  1  yq  rcL  }g  159  PVIEU 

10  PROMPT  S8  RCL  10  jig  x>V>  169  PSE 

11  STO  03  61  Xt2  111  QTO  "TPEa  161  RCL  18 

12  "PRS  SO.  LP""  62  ♦  112  i/=  162  VIEW  1? 

13  PROMPT  63  S9RT  1 13  R&.  89  163  PSE 

14  STO  04  6*  STO  11  114  *  164  ■RADIUS''" 

15  "PRS  EAST  LQTT  65  "GEO.  SPEED-"   115  qtuH  165  PROMPT 

16  PROMPT  66  BVIEH  U6  CHS  166  STO  25 

17  STO  05  67  PSE  « 17  270  167  8 

19  -PRS  HEST  LflT""  6P  RCL  11  113  +  i$8  EHTERt 

19  PROMPT  69  VIEW  II  119  §j(j  13  1^9  pcL  25 

29  STO  06  70  PSE  126  GTO  "GRQ?"  170  m=Y> 

2!  RCL  03  71  0  121*LBL  "TRE"  171  GTO  "PHTI ■ 

22  EHTERt  7?  ENTERS  122  1/X  172  RCL  11 
27  RCL  04  73  RCL  0?  133  pcL  09  173  RCl  25 

24  -  74  Y=V  ,04  *  {74  y 

25  100  ?5  GTO  "CHK"  125  hTwH  175  RCL  08 
li  *  76  0  126  Chs  176  / 

27  2.22  E05  7?  EHTERt  12?  270  177  .25 

23  /  ?3  RCL  10  128  +  178  «• 

29  1.3  79  X= V?  129  STO  IS  179  S8RT 

M  /  &   GT°  'cm'  130  GTO  "GRPD"  188  ^5 

31  CHS  31  0  131  ♦LBL  "CHIC'  131  * 

32  STO  07  32  EHTERt  132  fi  153  STO  2* 

33  RCL  01  83  RCL  09  133  EHTERt  133  Pu  1! 

34  EHTERt  84  X>Y?  134  RCL  Jg  154  ?Cl  2b 

35  SiH  35  GTO  "HSH"  135  y<=yv  195  , 
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42  /  92  RCL  09  J42  STO  18  192  RCL  08 

43  ST0  09  93  *  143  gto  "GRPD"  193  / 

44  RCL  05  94  mT&n  144*LBL  "CHKS"  194  STO  20 

45  EHTERt  95  STO  IS  145  a  195  RCL  20 

46  RCL  06  96  GTO  "GRAB"  146  EHTERt      196  ENTERt 

47  -  97*LBL  -TWO"  147  RCL  09      197  .25 

48  100  98  U't  J48  K<=Yv       198  + 

49  *  99  RCL  09  149  qtu  "EPST"    199  STO  28 

50  2.22  E05  100  *  150  270.0       200  0 


E .        Program 


Table  I 

Surface 

Winds  Program. 

i*LBL  "wIND" 

/ 

Jjjjl  pTUN 

151  STO  IS 

FIX  0 

1.3 

182  188 

152  GTO  ■GRPD" 

"CTft  LAT?" 

/ 

183  t 

153*LBL  "EPST- 

iia  KKOMPT 

54 

PCL  98 

184  STO  IS 

154  98.8 

95  STO  8! 

cc 

195  f-Tfi  T.pcn- 

155  STO  13 

96  "CTR  L0H':"" 

56 

STO  18 

186*LBL  "ySH" 

156  GTO  "GRPD" 

87  PROMPT 

c-> 

RCL  89 

197  8 

157*l8L  "GRPD" 

8?  STO  82 

58 

ENTERt 

188  ENTEPt 

158  "GEO.  MR  FROM 

8?  "PRS  NO.  LflT?" 

5* 

Xt2 

189  RCL  !8 

159  PVIEW 

ifl  PROHPT 

68 

RCL  18 

118  W> 

168  PSE 

11  STO  83 

61 

Kt2 

111  GTO  "TRE' 

161  PCL  IS 

12  "PRS  SO.  LPT'" 

62 

♦ 

112  1/X 

162  VIEW  18 

13  PROHFT 

63 

SSRT 

113  RCL  89 

163  PSE 

14  STO  84 

64 

STO  11 

114  * 

164  "RPDIIJS?" 

15  "PRS  EPST  LQT? 

■  65 

"GEO.  SPEEI 

:"   115  PTPN 

165  PROMPT 

16  PROMPT 

66 

PVIEH 

116  CHS 

166  STO  25 

17  STO  85 

67 

PSE 

117  278 

167  9 

18  "PRS  UEST  LOT? 

*  68 

PCL  11 

113  + 

168  ENTERt 

19  PROMPT 

69 

VIEM  11 

119  STO  18 

169  RCL  25 

28  STO  96 

78 

PSE 

128  GTO  "GRPIi* 

179  X<=V> 

21  RCL  83 

71 

H 

12HLBL  "IRE" 

171  GTO  "flNTI- 

22  ENTERt 

70 

ENTERt 

122  1/X 

172  PCL  11 

23  RCL  84 

i  0 

RCL  89 

123  PCL  89 

173  RCL  25 

24  - 

it 

W> 

124  * 

174  / 

25  188 

75 

GTO  "CHK" 

125  PTPN 

175  RCL  83 

26  * 

76 

9 

126  CHS 

176  / 

27  2.22  E85 

it 

<  i 

ENTERt 

127  278 

177  .25 

'??.  / 

78 

PCL  18 

123  t 

178  + 

?o  i  7, 

79 

X=V> 

129  STO  18 

179  SPRT 

3?  / 

88 

GTO  -CHKS" 

130  GTO  "GPfiD" 

188  .5 

31  CHS 

81 

i? 

131*LBL  "CHK" 

131  + 

32  STO  87 

32 

ENTERt 

132  8 

1S2  STO  26 

33  RCL  81 

83 

RCL  89 

133  ENTERt 

133  RCL  11 

34  ENTER? 

34 

r.)V 

134  RCL  18 

134  RCL  26 

35  -SIN 

85 

GTO  -USN" 

135  X<=V^ 

135  / 

36  1.4584  E-84 

36 

8 

136  GT0  "NOR" 

186  STO  27 

37  * 

37  ENTER? 

137  188.8 

137  GTO  "PCT" 

33  STO  88 

88 

RCL  18 

133  STO  18 

183*LSL  "PNTI" 

39  RCL  87 

39 

X>'i"' 

139  GTO  'GRPD* 

139  RCL  11 

48  ENTEPt 

98 

GTO  -TWO- 

14m*LBL  "NOR" 

198  RCL  25 

41  RCL  88 

91 

1/X 

141  368.8 

191  / 

42  / 

92 

PCL  89 

142  STO  13 

192  RCL  88 

43  STO  89 

93 

* 

143  GTO  "GRQD" 

193  / 

44  PCL  85 

94 

PTPN 

144»LBL  "CHKS" 

194  STO  29 

45  ENTERt 

95 

STO  IS 

145  9 

195  RCL  28 

46  RCL  86 

96  GTO  "GPPD" 

146  ENTERt 

196  ENTERt 

47  - 

97*LBL  *TUO- 

147  RCL  99 

197  .25 

48  188 

98 

1/X 

148  X<=Y-> 

198  + 

49  * 

99 

PCL  89 

149  GTO  "EAST" 

199  STO  28 

58  2.22  E95 

188 

* 

158  278.8 

288  8 

III.  Determining  the  Probability  of  Precipitation  for  Monterey,  California 

A.  Objective 

To  determine  the  probability  of  24-h  rainfall  for  a  regional  area 
(Monterey,  California)  for  the  period  0800  PST  forecast  day  to  0800  PST  of 
the  following  day,  for  the  wet-season  months,  November  through  April.  This 
procedure  is  based  on  Renard's  (1972)  honogram  method  and  statistics. 

B.  Principles 

General  Forecast  Procedure.  The  probability  of  24-h  rainfall  is 
determined  by  values  of  geostrophic  relative  vorticity  and  sea-level  pressures 
at  different  locations. 

A  first  forecast  parameter  is  obtained  by  taking  the  sea-level  pressure 
at  Eureka,  CA  (PEureka)  and  the  difference  between  the  500  mb 
geostrophic  relative  vorticity  at  Monterey  (  Cm)  and  a  point  8°  latitude 
upstream  from  Monterey  (  Cs).  Enter  these  values  into  a  pre-established 
statistical  table,  which  is  stored  on  the  calculator,  to  arrive  at  a  forecast 
variable,  Yi .   (See  Figure  4.) 

A  second  forecast  parameter,  Y2  ,  is  the  value  that  results  from  taking 
the  sea-level  pressure  difference  between  Monterey  and  Eureka  (APM-Eureka) 
and  Monterey  and  Las  Vegas  (APm-lv).  (See  Figure  5).  These  two  parameters, 
Yi  and  Y2 ,  are  entered  in  a  third  table  (Figure  6)  which  results  in  a  final 
value  that  corresponds  to  a  probability  of  average  rainfall  for  the  Monterey 
Peninsula. 

500  mb  Geostrophic  Relative  Vorticity.  To  calculate  the  500  mb 
geostrophic  relative  vorticity,  five  500  mb  heights  are  needed  to  use  finite 
differencing  technique  to  solve  (Cm  =  g/f  V2  Z)  where  m  is  the  geostrophic 
relative  vorticity  at  Monterey,  V2Z  is  the  Laplacian  of 'the  500  mb  height. 
For  example,   in  Figure  3  at  Monterey  (Zm  =  558  dm)  the   four  500  mb  heights 
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are  552  dm  (Zn :  500  mb  height  to  the  north),  566  dm  (Za ) ,  560  dm  (Ze),  and 
563  dm  (Zw)  with  a  grid  spacing  of  6°  latitude.  This  results  in  a  geostrophic 
relative  vorticity  of  1.8  x  10"5  sec-1. 


Figure  3 

Computation  of  geostrophic  relative 
vorticity  for  Monterey  (?  m)  and  a 
point  8°  latitude  upstream  (Cs). 


556    »t 
*>4     \ 

\VlaV 

Zv\        \ 

?fe\Zc 

*no  \ 

Zn^ 

v^ 

3m        73, 

1     ^^^ 

zs 

The  500  mb  geostrophic  relative  vorticity  is  computed  for  a  point  8° 
latitude    upstream  (8°  corresponds  to  a  typical  daily  progression  of  short- 
waves  in  the   westerlies  of  22  kt  during  the  winter  season)  from  Monterey  in 
the  same  way.    The  difference  between  the  two  vorticities  (£s°  -Cm)  is  used 
as  a  forecast  factor. 

Sea-Level  Pressure.  Yi ,  a  first  forecast  parameter,  is  obtained  from 
the  computed  vorticity  difference  and  a  sea-level  pressure  at  Eureka, 
California  (Figure  4).  Y2 ,  a  second  forecast  parameter,  is  obtained  from 
the  difference  between  the  sea-level  pressure  at  Monterey  minus  the  sea— level 
pressure  at  Eureka  (APM-Eureka )  and  the  sea-level  pressure  at  Monterey  minus 
the  sea-level  pressure  at  Las  Vegas  (APm-lv)  (Figure  5). 

Probability  of  Rainfall.  From  the  two  forecast  parameters,  Yi  and  Y2 ,  a 
value  Y3  is  obtained  which  corresponds  to  a  probability  of  rainfall  for  the 
Monterey,  CA  area  for  the  period  0800  PST  forecast  day  to  0800  PST  of  the 
following  day,  for  the  wet-season  months,  November  through  April.  (Figure  6.) 
The  corresponding  values  between  Y3  and  a  probability  of  average  rainfall  are 
on  data  cards  that  are  entered  during  the  program  execution. 


12 


,**-,<«*. 

9T, 

3  IO0O 

IO02. 

/o<y+ 

""Ofc 

loot 

o'OlO 

(Oil 

lOl^ 

/O/fcC'-'f 

10 10 

107  : 

27 

l?5 

7  5J  3.)      O-O  |oC    tOO   |00 

d-s 

2   1  |  4  i    **  <, 

j>'        1 
2?bo|  00 

0  0 

24 

90 
21 

I'M 

*         ! 

05 

\oo  |  0  0 

0.0 

1 

.'■7 

3  1 

3.1 

5.8 

«»•©  |  0  O 

0  0 

O  0 

1 

3..0 

\ 

i?5 

\ 
7  5 

3  2 

\ 

0  O 

00 

I 

f! 

b  O 

TO 

S  4     '  ■  & 

j 

>oo 

00 

18 

\ 
37.6 

C50 

\ 

5=1 

I- 7  |  os<x. 

e.or 

5  « 

7   V 

9  2 

.■7J3  1 

\ 
00 

0  a 

IS 

\ 

\ 

in   r\ 

4  u 

k    " 

/•" 

1  /. 

53 

1  u  ->, 

/ 

10  «* 

12    1 

HU7 
•   \ 

<->.H 

mY) 

'l1 

\ 
57  3  < 

JH>~7 

*7* 

\ 

7-5  J  SO 

4  3 

it? 

r 

14^3 

j     \ 
JS^P   14.1      i  v,' 

1  S 

\ 
CO 

1 

fcX> 

WV0 

,37.9 

r 

2J..3J  l\0 

- 
14 

7f- 

M'3 

> 

>71 

22.5 

24.4 

7  <• 

2  7 

1 
O  6J 

b 

wo 

(.2  5 

s 

,59  0 

27.5 

22"  6 

"2J-0 

~zl  6 

27  9 

'3/  O 

32. G»  22.0 

14 

S  4. 

1 

\ 
09 

3* 

■fcO-7 

(.33 

(>3.0 

'#0 

45^0 

5"?0 

37  5 

4/3 

434 

35  0 

3/  O. 

10  V 

'°/l 

4o 

Ota 

0 

%S 

52  5 

575 

<*3 

too 

M.2 

52. 3 

41*5 

v44.l 

3SO 

/ 

•270 

ita-O 
/ 

7  / 

4  2 

0  4 

-3 

JO 

-b 

41.3 

4^0 

■g-1 

S25 

•41.  3 

55.4 

ta<r  fc 

45,0' 

J 

'Hl-O 

/ 

3pO 

w.§ 

V^fo 

3  ? 

1-4- 

2SO 

Vo 

43.0 

350 

375 

37  9 

375 

32S 

27.  q 

\<t  $. 

'U.3 

/ 
fc.l 

21 

2.7 

07 

0 

"1 

12.% 

15^2 

11  0 

22  5 

2T0 

"l5~7 

T^.7 

2I-4 

17,1 

'11  3 

b/fc, 

27 

0  s 

0t>' 

DC 

-12 

b-+ 

C.7" 

iTV 

*i3~r 

I5-0- 

I5"D 

'l2  5 

9  0 

1.*f 

/■<? 

ot 

0  0 

OO 

00 

-IS 

45 

4-7 

S3 

(.0 

fc~4 

"ff-r 

7-5" 

fe-4 

3^2 

• 

• 

OO 

0  0 

OO 

00 

-l« 

21, 

2  S 

31 

3-6 

3  e 

3-4- 

2  9 

I.I 

erf 

00 

00 

OO 

0.0 

1 

00 

AP. 


M-LV 


-10          2.           4-           U          %          lO,  .,12-        1  + 

10 

% 

u 

e 
HO 

2 

O 

-2 
0 

-4 

-fc 

-t 

-» 

-1  z 

0 

-14 

0.0 

00 

0  0 

■o-7 

3  1 

£     7 

/ 

7ta 

I 
1  3    '0-1 

OO 

OO/ 

lb 

4-3 

1.  7 

1 

o,o- 

14 

50 

1< 

10-4 

12  '• 

.».« 

15 

UJ> 

ip'U 

if4 

(60^ 

.11 1 

';■/  j 

1" 

:2,0'|24;(i 

50/ 

II. 3 

i.ta.r 

"21.  u 

23-4 

i-C,    i 

0  4|:;3.«j4l.5 

'    1 

p.  1 

/ 

2;^ 

3l.V 

3H  t 

J"    1  ^'; 

r 

5A 

1 

'2tS 

31.0 

34.0 

iTJ'ii.arhs.yy-ui  h 

o«e 

— T — 
20  O 

22^(, 

3oN0 

n  '• 

•»»  5 

tc'o  I  l^f  CJ 

00 

oi>  3  « 

>* 

> 

-. 

o.« 

10.7^ 

n  0 

i\  0  'iVT_ 

."<,OJ.-5  7»tl  0 

o.p, 

1,7  "a 

"^ 

4i  r 

•/«!  5 

5~7    5 

3^'  ief5i  7*i 

*i^' 

<J7<5 

>/'o 

u.4 

ii  : 

^■J.n 

5r7f'+y.3 

., 1 .' 

.  1  / 

W-S 

W.I 

fc+  1 

U  5 

ta0.7U:,| 

lao.%so»'.WSt<1                                                             '•o 

AP, 


Eureka 


M-Eureka 


Figure  4 

Yi  =  ffPEureka,  Ca  -Cm):  first  forecast 
factor  of  the  raincaster. 
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Figure  5 

Y2  =  fCAPM-cur^ica  ,  APm-lvj: 
second  forecast  factor 
the  raincaster. 


Figure  6 

Yt    =  f(Yi  ,    Ys  >:    third 

forecast    factor  of  the 
raincaster. 
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c. 

1) 

2) 

3) 
4) 
5) 
6) 

7) 
8) 
9) 
10) 


User's  Guide 

Make  sure  the  printer  is  attached  and  in  the  NOR  mode.   Both  printer  and 

calculator  should  be  on  and  the  calculator's  memory  cleared. 


XEQ 


alpha   SIZE   alpha 


will  prompt  SIZE 

Type  035. 

Load  the  program  cards  (1-18).   Continue  after  the  prompt  "WORKING." 


XEQ 


alpha  |  RAIN   alpha 


Enter  500  mb  heights   (dm)  for  each  prompt.   (For  example,  ZN:  8°  means 
height  at  north  position  at  a  point  8°  latitude  upstream  from  Monterey. ) 


Hit  R/S  after  each  entry. 


Enter  pressure  values  for  Monterey,  Las  Vegas,  and  Eureka. 
Hit 


R/S  after  each  entry. 


11) 


A  data  card  prompt  will  occur  after  all  the  information  has  been  entered. 

For  the  first  time,   use  data  set  #1;  the  second  time,  use  data  set  #2, 

etc. 

The  final  number  that  appears  after  all  data  card  insertion  instructions 

is  the  probability  of  average  rainfall  for  the  Monterey  Peninsula. 

NOTE:   If  the  program  ends  with  "OUT  OF  RANGE",  then  the  values  used  are 

either  too  small  or  too  large.   Try  another  set  of  values  that 

will  fit  in  the  stored  data  sets. 
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D.   Examples 

Example  A  illustrates  the  use  of  the  program  "RAIN"  in  determining  the 
probability  of  precipitation  at  Monterey,  CA.  At  1200Z  25  November  1981  for  a 
point  8°  latitude  upstream  from  Monterey,  the  500  mb  heights  (in  decameters) 
6°  latitude  to  the  north,  east,  south,  west  and  at  the  center  point  are  541.0, 
543.0,  562.5,  563.0  and  555.0  respectively.  (See  Figure  7.)  Around  Monterey, 
the  500  mb  heights  are  550.0  dm  6°  latitude  to  the  north,  555.0  dm  at  the 
eastern  location,  573.0  dm  to  the  south,  576.0  dm  to  the  west,  and  567.0  dm  at 
Monterey.  These  heights  are  used  to  obtain  the  500  mb  relative  geostrophic 
vorticity  values  used  in  generating  forecast  variables.  The  1200Z  sea-level 
pressures  at  Monterey,  Las  Vegas  and  Eureka  (1018.8  mb,  1006.2  mb,  1017.3  mb) 
are  also  required  to  compute  forecast  variables.  (See  Figure  8a.)  These 
forecast  variables  result  in  a  probability  of  precipitation  of  90%.  Rain  did 
verify  during  the  24-h  period,  as  noted  from  the  1200Z  26  November  1981 
observation  (Figure  8b.) 
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Figure  7 
500  mb  analysis  for 
1200Z  25  November  1981 
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Figure  8a 
Surface  observations 
for  1200Z  25  November  1981 


Figure  8b 
Surface  observations 
for   1200Z  26  November  1981 


Table  II 
Example  of  the  Probability  of  Precipitation  for  Monterey,   California 
Program. 


ENTER  ZN:S 


ENTER  ZES 


ENTER  ZS:8 


ENTER  ZU:S 


ENTER  Z*:S 


ENTER  ZM:H 


ENTER  ZE:fl 


ENTER  ZS:N 


ENTER  ZW:H 


ENTER  Z*:M 


54 1. 


5*3. 


562.5 


563. 


555. 


548. 


554. 


569.5 


579.5 


ENTER  PRE3  HOMT 

1*818.8 
ENTER  PRES  LVGS 
I »086. 2893 

ENTER  PRES  ERKfi 
1..  91 7.  3999 

23=  L018.0888 

X2=-33.9999 

X3=  2,9909 

X4=  14.9999 

INSERT  CRD  13 

Yl=-28.0088 
INSERT  CRD  3 

Y2=  25.9996 
90.9999 

Probability  of 
Precipitation  is  90%. 
Data  verified  from 
observed  case  at 
1200Z  26  November  1981 
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E. 


Program 


Table  III 
Probability  of  24-h  Precipitation  for  Monterey,  California  Program. 


8!*LBL  "RPIW- 

92  "ENTER  ZH8 

93  Pv-ORPT 
34  STO  98 

85  "ENTER  ZE:?' 
96  PROMPT 
87  RCL  88 
9?  + 

89  STO  88 
18  "ENTER  ZS' Pi- 
ll PROMPT 
12  RCL  88 
.13  + 

14  STO  88 

15  -ENTER  ZW:y- 

16  PROMPT 
1?  RCL  89 

18  + 

19  STO  88 

28  "ENTER  Z8:8" 

21  PROMPT 

22  4 

23  * 

24  CHS 

25  RCL  88 

26  ♦ 

27  STO  m 

n   "ENTER  ZN:H' 

29  PROMPT 

38  STO  81 

31  "ENTER  ZE:M- 

32  PROMPT 
37  RCL  81 

34  + 

35  STO  81 

36  -ENTER  ZS:«" 

37  PROMPT 

33  RCL  81 

39  + 

48  STO  91 

41  -ENTER  ZH:M" 

42  PROMPT 

43  RCL  81 

44  + 

45  STO  91 

46  "ENTER  Z8M' 

47  PROMPT 

48  4 

49  * 

58  CHS 


51  RCL  8 1 

52  + 
57  CHS 

54  RCL  8? 

55  + 

56  STO  82 

57  -ENTER  PRES  MONT' 

58  PROMPT 

59  XEQ  "ROUND" 
68  STO  83 

61  "ENTER  PRES  LVGS" 

62  PROMPT 

63  XEQ  -ROUND" 

64  STO  84 

65  'ENTER  PRES  ERKfl- 

66  PROMPT 

67  XEQ  "ROUND" 

68  STO  81 

69  RCL  81 
79  2 

71  MOD 

72  X=8^ 

73  GTO  35 

74  RCL  81 

75  1 

76  + 

77  STO  81 
73  GTO  35 
79HBL  85 
38  RCL  81 

31  "Xl=" 

32  PCP 

33  PCX 

34  PDV 

35  RCL  82 

36  XEQ  -ROUND" 
87  STO  82 
88HBL  77  . 

89  PCL  82 

98  3 

91  MOD 

92  X=8? 

93  GTO  86 

94  1 

95  X=Y? 

96  GTO  83 

97  GTO  87 
98*LBL  S3 

99  RCL  82 
188  1 


181  - 

182  STO  82 
193  GTO  77 
184*LBL  37 

185  RCL  82 

186  1 
197  + 

188  STO  82 

189  GTO  77 
119H6L  86 

111  RCL  82 

112  -X2=- 

113  PCS 

114  PCX 

115  PDV 

116  RCL  93 

117  RCL  81 

118  - 

119  XEQ  -FIXEVEN- 
129  RCL  85 

121  STO  88 

122  "X3=- 
!23  PCP 

124  PCX 

125  PDV 

126  PCL  83 

127  RCL  84 

128  - 

129  XEQ  "FJVEVEN- 
138  RCL  85 

131  STO  84 

132  "X4=- 

133  PCP 

134  PCX 

135  PDV 

176  994 

137  STO  87 

138  27 

139  STO  83 
149  2 

141  STO  89 

142  RCL  89 

143  STO  83 

144  RCL  91 

145  GTO  48 
146H8L  48 
147  RCL  87 
143  X>Y? 
149  GTO  23 
159  X=Y? 


151  GTO  88 

152  XEQ  98 

153  X=Y? 

154  GTO  89 

155  XEQ  88 

156  X=Y? 

157  GTO  19 

158  XEP  88 

159  X=Y? 
168  GTO  11 

161  XEQ  98 

162  X=Yv 

163  GTO  12 

164  XEQ  88 

165  X=Y? 

166  GTO  13 

167  XEQ  88 

168  X=Y/ 

169  GTO  14 
173  XEQ  88 

171  X=Y? 

172  GTO  15 

173  XEQ  88 

174  X=Y? 

175  GTO  16 

176  FS^  89 

177  GTO  23 

178  XEQ  88 

179  X=Y? 
138  GTO  17 
131  XEQ  98' 

182  X=Y? 
1S3  GTO  13 
134  XEQ  88 
185  X=yv 

136  GTO  19 

137  XEQ  88 

183  X=Y? 
189  GTO  28 

198  XEQ  88 

191  X=Y? 

192  GTO  21 

193  XEQ  88 

194  X=V? 

195  GTO  22 

196  GT0  23 
197*LBL  88 
193  RCL  99 

199  + 
298  RTN 


291*LBL  88 

292  29.935 

293  "INSERT  CRD  !" 

294  BVIEH 

295  RDTPX 

296  FS"'  00 


.-  GTQ  V, 

288  FS?  81 

•700  f-7 


28?  GTO  44 

218  GTO  24 
21ULBL  8? 

212  28.835 

213  "INSERT  CRD  2" 

214  pv?Ffc 

215  RDTPX 

216  FS?  88 

217  GTO  33 
213  FS?  91 

219  GTO  44 
228  GTO  24 
221*LBL  10 

222  20.035 

223  "INSERT  CRD  3" 

224  PVIEH 

225  RDTPX 

226  FS?  09 

i-i-i     b'.'J    i-j 

228  PS'}  01 

229  GTO  44 
238  GTO  24 

"3 1  i  A I  CI       *  * 

CSlHoL    il 

232  28.835 

233  "INSERT  CRD  4- 

234  PVIEU 

235  RDTPX 

236  FS"'  88 

237  GTO  33 

23g  f§v  pi 

23?  GTO  44 

248  GTO  24 
24l*LBL  12 

242  28.835 

243  "INSERT  CRD  5" 

244  PVIEW 

245  RDTPX 

246  XROM  23.33 

247  RCL  81 
243  FS?  8i 

249  QTO  44 

250  GTO  24 
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251*LBL  13 

252  28.835 

253  "INSERT  CRD  6* 

254  flVIEH 

255  RBTflX 

256  flVIEH 

257  XPOH  28.33 

258  GTO  33 

259  FS?  81 

260  GTO  44 

261  GTn  24 
262»LBL  14 

263  29. 835 

264  -INSERT  CRD  ?' 

265  AVI  EH 

266  RBTflX 

267  XRGH  28,33 

268  GTO  33 

269  FS?  91 
278  GTO  «4 
271  GTO  24 
272*LBL  15 

273  28.835 

274  "INSERT  CRD  8" 

275  SVIEH 

276  RDTflX 

277  GTO  33 

278  FS''  33   ■ 
279*LBL  88 
288  GTO  44 
281  GTO  24 
232*LBL  16 

283  28.835 

284  'INSERT  CRD  9" 

285  PVIEH 

286  RDTflX 

287  XROH  28.33 

288  RCL  81 

289  FS?  81 

298  GTO  44 
291  GTO  24 
292HBL  17 

293  28.835 

294  -INSERT  CRD  18" 

295  AVI  EH 

296  RDTfiX 

297  mm   23,33 
298*LBL  88 

299  GTO  44 
388  GTO  24 


38ULBL  13 

382  28.835 

383  "INSERT  CRD  11" 

384  flVIEH 

385  RDTAX 

386  XROH  28,33 

387  X>8':' 

388  RCL  12 

389  GTO  24 
318*L8L  19 

311  28.835 

312  -INSERT  CRD  12" 

313  AVIEH 
3! 4  RDTAX 

315  GTO  44 

316  XROH  28.33 

317  GTO  24 
318*LBL  28 
319  28.835 

328  "INSERT  CRD  13" 

321  flVIEH 

322  RDTAX 

323  FS?  81 

324  GTO  44 

325  GTO  24 
326*LBL  21 
327  28.835 

323  "INSERT  CRD  M" 

329  AVI  EN 

338  RDTAX 

331  flVIEH 

332  XROH  28.33 

333  GTO  44 

334  GTO  24 
335*LBL  22 

336  28.835 

337  "INSERT  CRD  15* 
333  flVIEH 

339  RDTAX 

348  GTO  d 

341  XROM  28,33 

342  RCL  12 

343  GTO  24 
344*L8L  24 

345  RCL  82 

346  RCL  88 

347  X=Y? 
343  GTO  58 

349  XEQ  25 
358  X=Y? 


351  GTO  51 

352  XEQ  25 

353  x=V? 

354  GTO  52 

355  XEQ  25 

356  X=Y? 

357  GTO  53 

358  XEQ  25 

359  X=Y? 
368  GTO  5* 

361  XEQ  25 

362  X=Y? 

363  GTO  55 

364  XEQ  25 

365  X=Y? 

366  GTO  56 

367  XEQ  25 

368  X=Y? 
363  GTO  57 
378  XEQ  25 

371  X=Y? 

372  GTO  58 

373  XEQ  25 

374  X=Y? 

375  GTO  59 

376  XEQ  25 

377  X=Y? 

378  GTO  68 

379  m  25 
388  X=^ 
331  GTO  61 

382  XEQ  25 

383  X=Y? 

334  GTO  62 

335  XEQ  25 

336  X=Y? 
387  GTO  63 
333  XEQ  25 
339  X=Y? 
398  GTO  64 

391  XEQ  25 

392  X=Y? 

393  GTO  65 

394  GTO  26 
395*LBL  25 

396  3 

397  - 

398  RTN 
399*LBL  26 

488  XEQ  -ROUND" 


481  STO  86 

482  5 

483  HOD 

484  STO  13 

485  X=mv 

486  GTO  31 

487  RCL  13 

488  1 

489  X=r'> 
418  GTO  38 

411  RCL  13 

412  2 

413  X=V? 

414  GTO  29 

415  RCL  13 

416  3 

417  X=Y? 

418  GTO  28 

419  RCL  13 
428  4 

421  X=Y? 

422  GTO  2? 
423*LBL  27 
424 
425 
426 


RUL 
1 


Hh 


j-*i" 

*♦-• 


96 


+ 
CTfl 

423  GTO  31 
423*LBL  28 
438  RCL  86 

431  2 

432  + 

433  STO  86 

434  GTO  31 
435*LBL  29 

436  RCL  86 

437  2 

438  - 

439  STO  86 
448  GTO  31 
441 ♦LBL  38 

442  RCL  86 

443  1 

444  - 

445  STO  86 

446  GTO  31 
447*LBL  31 

448  FS'  88 

449  GTO  32 
458  FS?  81 


451  GTO  35 

452  RCL  86 

453  STO  18 

454  "¥{=" 

455  PCS 

456  AC* 

457  flDV 

45P  -? 


459 


!!  H, 


468  18 

461  STO  88 

462  2 

463  STO  89 

464  SF  88 

465  RCL  83 

466  GTO  4h 
467*lBL  32 

468  RCL  % 

469  STO  12 

478  "Y2=" 

471  RCfi 

472  flCX 

473  flDV 

474  8 

475  STO  87 

476  78 

477  sto  88 

473  5 

479  STO  89 

438  SF  81 
481  CF  88 

432  RCL  18 

433  GTO  48 
434*LBl  23 

435  FS?  8! 

436  GTO  -END" 

437  "OUT  OF  RANGE" 
433  QVIEH 

439  RTN 
498*LBl  "END" 

491  FS?  82 

492  GTO  99 

493  SF  82 

494  RCL  86 

JijC  ', 

n  ?-.'  1 

496  ST0  88 

497  XEQ  88 
493  X>Y? 

499  GTO  66 

598  XEQ  88 
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591  X>Y? 

551 

48 

681  X=v? 

592  GTO  67 

552 

X>Y? 

692  GTO  61 

583  XEQ  88 

553 

GTO  83 

693  XEQ  88 

584  X>Y? 

554 

GTO  84 

£,M4  X=Y^ 

585  GTO  68 

555 

RT* 

695  GTO  62 

596  XEQ  88 

J  JO 

►LSI  99 

696  GTO  63 

58?  X>Y? 

re-* 
J-.U 

PCX 

687  GTO  26 

588  GTO  69 

j  36 

PBV 

698*LBL  35 

589  XEQ  89 

559 

CF  81 

699  RCL  86 

518  JOY? 

568 

CF  82 

619  "v3=" 

511  GTO  78 

561 

RTN 

611  PCP 

512  XEQ  98 

562' 

►LBL  33 

612  PCX 

513  X>Y? 

563 

_■■» 

613  PBV 

514  GTO  71 

564 

STO  89 

614  GTO  23 

515  XEQ  88 

565 

RCL  94 

615*>LBL  44 

516  X>Y? 

566 

RCL  88 

616  -5 

517  GTO  72 

56? 

a-  i ' 

617  STO  99 

518  XEQ  88 

568 

GTO  58 

613  RCL  12 

519  X>V> 

569 

XEQ  99 

619  RCL  88 

528  GTO  73 

578 

X=T'' 

629  X=Y? 

52!  XEQ  88 

571 

GTO  51 

621  GTO  59 

522  X>Y? 

J  t  C 

XEQ  98 

622  XEQ  88 

523  GTO  74 

573 

X=Y? 

623  X=Y? 

524  XEQ  89 

574 

GTO  52 

624  GTO  51 

525  X>Y? 

575 

XEQ  88 

625  XEQ  89 

526  GTO  75 

576 

x=v? 

626  X=Yv 

527  XEQ  88 

J  !  1 

GTO  53 

627  GTO  52 

528  X>Y? 

578 

XEQ  99 

628  '^£9   99 

529  GTO  76 

579 

X=Y? 

629  X=Y? 

538  XEQ  89 

588 

GTO  54 

638  GTO  53 

531  XEQ  88 

581 

XEQ  88 

631  XEQ  98 

532  X>Y^ 

532 

632  X=Y? 

533  GTO  78 

533 

GTO  55 

633  GTO  54 

534  'fi.iQ   99 

534 

XEQ  88 

634  XEQ  88 

535  XEQ  88 

585 

X=Y? 

635  X=Y? 

536  X>Y? 

536 

RCL  26 

636  GTO  55 

537  GTO  79 

537 

GTO  56 

637  XEQ  88 

538  XEQ  88 

538 

XEQ  88 

633  X=Y? 

539  XEQ  89 

539 

X=Y? 

639.  GTO  56 

548  X;>Y? 

598  GTO  5? 

648  XEQ  89 

541  GTO  38 

591 

XEQ  88 

64!  X=Y? 

542  '<£Q  88 

592  X=Y? 

642  GTO  5? 

543  XEQ  88 

593 

GTO  58 

643  XEQ  88 

544  X>Y? 

594 

XEQ  89 

644  X=Y? 

545  GTO  81 

595 

A-T  '. 

645  GTO  58 

546  CLX 

596 

GTO  59 

6^6  XEQ  88 

547  29 

597 

XEQ  99 

64?  X=Y? 

548  x>y-? 

598 

X=Y? 

643  GTO  59 

54?  GTO  82 

59? 

GTO  69 

649  XEQ  98 

558  CLX 

688 

XEQ  88 

659  W> 

651  GTO  68 

652  XEQ  98 

653  X=v" 

654  GTO  61 

655  XEQ  99 

656  X=V? 
65?  GTO  62 

658  XEQ  89 

659  K=Y? 
668  GTO  63 
66 i  XEQ  89 

662  X=Y? 

663  GTO  64 

664  GTO  65 

665  GTO  26 
666*LBL  59 
66?  RCL  29 
668  GTO  26 
669 ♦LBL  51 

678  RCL  21 
671  GTO  26 
672*LBL  52 

673  RCL  22 

674  GTO  26 
675>LBL  53 

676  RCL  23 

677  GTO  26 
678HBL  54 

679  RCL  24 
638  GTO  26 
681*L8L  55 
682  RCL  25 
633  GTO  26 
684+LBL  56 
685  RCL  26 
636  GTO  26 
687*LBL  57 

638  RCL  2? 

639  GTO  26 
698*L8L  53 

691  RCL  23 

692  GTO  26 
693*LBL  59 

694  RCL  29 

695  GTO  26 
696*LBL  68 
697  RCL  38 

693  GTO  26 
699+L8L  61 
788  RCL  31 


78!  GTO  26 
782*LBL  62 
783  RCL  32 
794  GTO  26 
?85*LBL  63 
796  RCL  33 
787  GTO  26 
788+LBL  64 
789  RCL  34 

718  GTO  26 
71! ♦LBL  65 

712  RCL  35 

713  GTO  26 
714*i_BL  66 

715  18 

716  GTO  23 
717*LBL  67 

i  1 0  1  J 

719  GTO  23 
728+LBL  63 

721  25 

722  GTO  23 
723»LBL  69 

724  38 

725  GTO  23 
726*L8L 


■fl 


728  GTO  23 
?29*LBL  71 

738  48 

73!  GTO  23 
732HBL  72 

733  45 

734  GTO  23 
735+L6L  ?3 
736  58 
73?  GTO  23 
738*LBL  74 

739  55 
748  GTO  23 
74HLBL  75 

742  68 

743  GTO  23 
744*LBL  76 

745  65 

746  GTO  23 
74?*LBL  78 

748  78 

749  GTO 
758*LBL 


^5!  75 
752  GT0  23 
753*lBL  88 
754  38 

f  j  J   'j  t  V  ■-  - 

756*lBl  31 
75?  85 
758  GTO  23 
759«>LBL  82 

768  99 

761  GTO  23 
762*LBL  33 

763  95 

764  GTO  23 
765*LBL  64 
766  188 
76?  GTO  23 
763*LBL  "ROUND" 

769  FIX  8 

778  RND 

771  FIX  4 

772  RTN 
7?3*LBL  -FIXEVEH- 

774  STO  85 

775  2 

776  MOD 
77?  X=8? 

773  RTN 

779  RCL  85 
739  1 

731  + 

732  STO  95 

733  RTH 

734  .END. 


L.-J 

7<J 
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IV.   Contrail  Formation 

A.  Objective 
To  determine  whether  or  not  contrails  will  form  at  a  certain  level  in  the 

atmosphere  based  on  the  curves  from  the  Contrail  Forecasting  Manual  (Chief  of 
Naval  Operations,  1964). 

B.  Principles 
The   critical  temperature  determines  if  contrails  will  form  at  a  level. 

A  least  squares  approximation  is  given  for  the  critical  temperature: 
Tent  =  ai  +  02    In  p  +  a3  (In  p)2  +  a4  RH  +  as  (RH)2; 
where : 

p  is  the  pressure  level; 

RH  is  the  relative  humidity; 

ai  =  -90.4994;  a2  =  3.4232;  aa  =  0.5587;  a4  =  -0.0372;  as  =  0.0012. 

This  equation  was  determined  from  the  curves  in  Figure  9  of  the  Contrails 
Forecasting  Manual   (Chief  of  Naval  Operations,  1964). 

If  the  ambient  temperature  is  less  than  the  critical  temperature,  then 
contrails  will  form.  If  the  ambient  temperature  is  greater  than  the  critical 
temperature,  then  contrails  will  not  form. 

There  is  a  narrow  range  of  temperatures  where  the  relative  humidity  is 
important.  If  the  dew  point  is  known,  the  relative  humidity  can  be  calculated 
by  using  the  definitions  of , saturation  vapor  pressure  (es)  and  relative 
humidity  (R.H.). 

es(T)  =  6.11  exp  (L/Rw  (1/273  -  1/T)) 

R.H.  =  100  *  (e/es) 

where  es (T)  is  the  saturation  vapor  pressure 
e  is  the  actual  vapor  pressure 
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T  is  the  temperature  (°K) 

L  is  the  latent  heat  of  evaporation  (J/kg) 

Rw  is  the  specific  gas  constant  for  pure  water  vapor  (J/kg  °K) 

If  the  dew-point  temperature  is  unknown,  the  relative  humidity  is 
estimated  and  a  +2°C  error  margin  is  assumed  when  computing  the  critical 
temperature.  If  the  ambient  temperature  is  within  +2°C  of  the  critical 
temperature,  a  "probably"  will  precede  the  "contrails"  or  "no  contrails" 
message.   The   following  criteria  are  used  to  estimate  the  relative  humidity: 

Tamb  <  Tent  — >  (RH  =  0)  — >  contrails 
Tamb  <  Tent  ~ >  (RH  =  100)  — >  no  contrails 
p  <  225  — >  (stratosphere)  — >  RH  -  0% 
p  >  300  — >  (not  in  upper  troposphere)  — >  RH  -  40% 
225  <  p  <  300  — >  (upper  troposphere)  — >  RH  -  40%    UNLESS 
cirrus  are  at  this  level  — >  RH  —  60% 
flow  is  known  and  from 

moist  region  — >  RH  -  60% 
dry  region  — >  RH  -  0% 
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c. 

1) 

2) 
3) 
4) 
5) 
6) 
7) 

8) 
9) 
10) 
11) 

12) 


User's  Guide 

Make  sure  both  the  calculator  and  printer  are  on  and  that  the  printer  is 
set  to  the  MAN  mode. 
SIZE 


XEQ   alpha 


alpha 


Prompts  SIZE ;  Type  100. 

Load  all  three  cards,  sides  1-6. 

After  "WORKING"  signal,  program  is  ready  for  execution. 
TRAILS 


XEQ 


alpha 


alpha 


Program  asks  about  CHANGING  FLAGS.   For  the  first  run  through,  leave  the 
flags  as  they  are  and  continue  by  hitting  the  \   R/S  key. 
Enter  a  pressure  level  in  mb. 
Press 


R/S  key.   (Hit  R/S  after  each  entry.) 


13) 


14) 


Enter  a  temperature  in  °C. 

A   prompt  for  a  dew— point  temperature  will  be  displayed.   If  Td   is 

known,  type  a  1;  if  Td  is  unknown,  type  a  0. 

The  program  will  then  decide  whether  or  not  contrails  can  be  formed  at 

that  particular  level  depending  upon  the  criteria  discussed  on  the 

preceding  pages.   If  the  dew-point  temperature  is  unknown,  there  might  be 

more  "yes"  or  "no"  questions  depending  upon  the  pressure  level  and 

ambient  temperature.  Answer  "1"  for  "yes"  and  "0"  for  "no." 

The  program  will  determine  whether  or  not  contrails  will  form  at  a 

certain  pre-established  pressure  level. 

After  the  answer  is  displayed,  execution  can  continue  for  another  level 


by  pressing  1  R/S| 
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By  changing  the  flags,  two  options  are  possible.  Setting  Flag  1  will 
make  the  prompts  shorter  and  less  specific.  Setting  Flag  2  means  that 
the  program  assumes  that  the  dew-point  temperature  is  known  and  skips  the 
questions  concerning  relative  humidity  estimation. 


To  change  a  flag:   a)   Press  gold  button  and  | SF  |  key. 

b)  Program  prompts  SF  . 

c)  Enter  either  (or  both):   01  -  short  prompt 

02  -  known  Td 


d)  To  clear  flags,  press  gold  button  and  [  CF  key. 

e)  Program  prompts  CF  . 

f)  Enter  01  or  02  to  clear  any  flag  that  was  set. 
The  program  assumes  all  flags  are  clear  each  time  it  is  first  executed. 

D.   Examples 

Examples  A-E  in  Table  IV  illustrate  the  use  of  the  program  "TRAILS"  in 
deciding  whether  contrails  will  form  from  upper-level  sounding  data.  In 
example  A,  where  no  flags  have  been  set,  a  pressure  of  450  mb  and 
temperature  of  -35° C  are  entered  in  the  "TRAILS"  program.  The  dew-point  is 
unknown  in  this  case  so  a  "0"  is  entered.  The  program  decides  that  contrails 
will  not  be  possible  with  these  conditions  which  agrees  with  the  contrail 
curve  diagram.  (See  Figure  9.)  Example  B  is  similar  to  A  (no  flags  set  and 
dew  point  unknown),  but  contrails  will  probably  form  with  the  particular 
upper-level  sounding  data.  In  example  C,  no  flags  are  set,  but  the  dew  point 
is  known.  With  a  pressure  of  250  mb,  a  temperature  of  -60.0°C,  and  dew  point 
of  -62.0°C,  contrails  will  form.  Examples  D  and  E  have  flags  set,  which 
result  in  shortened  prompts  for  both  cases,  and  for  example  E,  the  dew  point 
is  assumed  to  be  known.  Both  examples  D  and  E  illustrate  that  contrails 
will  not  form  in  these  environmental  conditions. 
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Table  IV 
Examples  of  the  Contrail  Formation  Program. 


P=?  MB 


458.8888 


-35,8808 
DO  YGU  KNOW 
BEHPOINT*  TD? 
YES=0NE 
K0=ZERC 
KNOW  TD? 

8.8888 
NO  CONTRAILS 

No  flags  set 
T,  unknown 


*** 


*** 


#** 


B 


P=?  MR 

258. 8888        *** 
T=?  C 

-54,88mm        *** 
DO  YQU  KNOW 
DEWPOIHT,  TD? 
YES=0NE 
N0=ZER0 
KNOW  W 

m.8888        *** 

ARE  there 

CIRRUS  AT 
THIS  LEVELS 
YES=0HE 
H0=ZER0 
CIRRUS? 

1.8888        *** 
PP08A8LY 
CONTRfilLS 

No  flags  set,  Td  unknown 


p=?  m 

225.8888 
fsO  f: 

-68.8888 
DO  Y9U  KNOW 
BEHPGINL  TP 
YES=0NE 
N0=ZERQ 
KNOW  TD? 

1.8888 
TD=?  C 

-62.8888 
CONTRAILS 


*** 


*** 


*** 


*** 


P*?  MB 

388.8888 
T=?  C 

-48.8888 
KNOW  TD? 

8.8888 
CIRRUS? 

8.8888 
KNOH  FLOW? 

1.8888 
MST  OR  DRY? 

1.8888 
HO  CONTRAILS 


#** 


*** 


*** 


*** 


*** 


*** 


No  flags  set 
T,  known 


Flag  01  set 
Unknown  T  A 


p=?  m 

858.8888 

T=?  C 

-18.8888 

TD=?  C 

-28.8888 

NO  CONTRAILS 

Flag  02  set 

T 

.  known 

*** 


*** 


*** 
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Figure  9 

Illustration  of  contrail  decisions:  curves  taken  from  the  Chief  of 
Naval  Operations  (1964). 
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E.       Program 


Table  V       Contrail   Formation  Program. 


91*LBL  " TRAILS ■ 
32  C'F  9  J 

83  CF  82 

84  • CHANGE  FLfiC 

85  PROMPT 
86*L6L  18 

87  TONE  1 

88  "P='  HB" 

89  FS'  2! 

18  PRfi 

11  PROMPT 

12  FS'  21 

13  PRX 

14  STO  88 

15  LN 

16  STO  81 

17  TONE  2 
15  "T=?  C" 

19  FS?  21 
'28  PRfi 

21  PROMPT 

22  FS'  21 

23  PRX 

24  STO  84 

25  FC?  82 

26  GTO  83 
27H8L  81 
28  XEQ  "RH" 
2?  XEQ  "TC- 
38*L8L  82 

31  RCL  84 

32  RCL  83 

33  X>Y? 

34  GTO  C 
35*LBL  D 

36  "NO  CONTRAILS 

37  TONE  5 

38  AVIEH 

39  ADV 
48  STOP 
41  GTO  18 
42*L6L  C 

43  "CONTRAILS" 

44  TONE  8 

45  fiVIEH 

46  ADV 

47  STOP 

48  GTO  18 
49*L8L  83 
58  FS?  81 


C"V 


51  GTO  U 

52  "DO  YGU  KNOW1 

53  AVIEH 

54  PSE 


55  -DEHPOIHT,  W 

56  fiVIEH 

57  PSE 

58  "YES=ONE" 

59  AVIEH 
68  PSE 

61  -NO=ZEPO" 

62  AVI  EH 

63  PSE 
64>LBL  11 

65  "KNOW  IB?" 

66  FS'  21 

67  PRA 

68  XROM  23.. 21 

69  XROH  28,33 
78  CHS 

71  1 

72  X=Y? 

73  GTO  81 

74  8 

75  STO  82 

76  XEQ  -TC" 

77  RCL  84 

78  RCL  83 

79  X>Y? 

38  GTO  C 
8i  188 

32  STO  82 

83  XEQ  'TC- 

84  RCL  83 

85  RCL  84 
36  X>Y? 

87  GTO  H 

88  8 

39  STO  82 
98  RCL  88 

91  225 

92  X>Y? 

93  GTO  A 

94  48 

95  STO  82 

96  388 

97  RCL  88 

98  X>Y? 

99  GTO  A 
188  68 


181  STO  82 

182  FS?  81 

183  GTO  12 

184  -ARE  THERE" 

185  AVIEH 

186  PSE 

187  -CIRRUS  AT- 

188  AVI  EH 

189  PSE 

118  -THIS  LEVEL0 

111  AVIEH 

112  PSE 

113  "YES=ONE" 

114  AVIEH 

115  PSE 

116  "NO=ZERO* 

117  AVIEH 

118  PSE 
H9*LBL  12 
128  "CIRRUS'- 

121  FS?  21 

122  PRA 

123  FS?  21 

124  XROM  28,33 

125  PRX 

126  X=Y? 

127  GTO  A 

128  XROM  23.33 

129  48 

138  STO  82 

131  FS'  81 

132  GTO  13 

133  -DO  YOU" 

134  AVIEH 

135  PSE 

136  -KNOW  FLOH?" 

137  AVIEH 

138  PSE 

139  "YES=GNE" 
148  AVIEH 

141  PSE 

142  "NO=ZERO" 

143  AVIEH 

144  PSE 
145*LBL  13 

146  -KNOH  FLOH?" 

147  FS?  21 

148  PRA 

149  PS'  21 
158  PRX 


151  XROM  28.33 

152  X=0? 

153  GTO  A 

154  FS?  81 

155  GTO  14 

156  "PROM  HOIST" 

157  AVIEH 
153  PSE 

159  "REGION- 

168  AVIEH 

161  PSE 

162  -=  ONE" 

163  AVIEH 

164  PSE 

165  -FROM  DRY" 

166  AVIEH 

167  PSE 

163  "=  ZERO" 

169  AVIEH 
178  PSE 
17ULBL  14 

172  "MST  OR  DRY'" 

173  FS?  21 

174  PRA 

175  XROM  28,33 

176  5 

177  PRX 

178  * 

179  GTO  39 
138  RCL  81 
131  X<8' 
182*L6L  P 

183  XEQ  "TC- 

184  RCL  83 

185  RCL  84 
136  - 

187  ABS 
133  2 
189  X<=Y? 
198  GTO  82 

191  "PROBABLY- 

192  AVIEH 

193  PSE 

194  GTO  82 
195*lBL  "RH" 

196  RCL  84 

197  273 

198  + 

199  XEQ  "ES" 
288  ST0  62 


201  TCNE  3 

202  "TD=?  C" 
283  FS'  21 
294  PRA 

285  FS?  2! 

286  PRX 

287  -XROM  28,: 

288  273 

289  * 

218  XEQ  "ES" 

211  RCL  82 

212  / 

213  188 

214  * 

215  ST0  02 

216  RTN 
217H8L  "ES" 

218  1/X 

219  CHS 

228  273 
221  1/X 

000  f 

223  2568 

224  * 

225  0.461 

226  / 

227  EtX 
oop  c 

229  * 

238  RTN 

231 ♦LBL  "TC" 

232  -90.4994 

233  ST0  03 

234  .5587 

235  RCL  81 

236  * 

237  3.4232 

07;:;  + 

239  RCL  81 

240  * 

241  ST+  03 

242  .0812 

243  RCL  82 

244  * 

245  -.8372 

246  ♦ 

247  RCL  02 

248  * 

249  ST+  03 
258  RTN 
251  .END, 
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V.   Determination  of  Lifting  Condensation  Level 

A.  Objective 

Given  a  surface  temperature,  surface  pressure  and  specific  humidity,  this 
program  calculates  the  lifting  condensation  level. 

B.  Principles 

The  lifting  condensation  level  (LCL)  is  the  level  to  which  a  parcel  of 
air  can  be  lifted  dry  adiabatically  before  it  becomes  saturated.  During  the 
lifting  process,  the  potential  temperature  of  the  air  parcel  and  the 
saturation  mixing  ratio  remain  constant.  The  actual  mixing  ratio  decreases 
and  eventually  equals  the  saturation  mixing  ratio.  The  lifting  condensation 
level  is  obtained  when  the  actual  vapor  pressure  equals  the  saturation  vapor 
pressure  (e  =  es ) .  Integrating  the  Clausius-Clapeyron  equation  from  273  K  to 
air  temperature  (T) 


T         T 

f  det  f     L  dT 

-L  es      }    Rw   T2 


where  es  is  the  saturation  vapor  pressure;  L  is  the  latent  heat  of 
condensation  (2.5  x  106  J/kg) ;  Rw  is  the  gas  constant  for  water  vapor  (462 
J/kg  °K);  and  T  is  the  temperature  (°K),  results  in  the  following  expression 
for  the  saturation  vapor  pressure: 

es  =6.11  exp  (- —  (---  -  -)). 
Rw   273   T 

Combining  the  hydrostatic  equation  and  the  ideal  gas  law  and  expressing  the 

temperature  as  To  -  Y  z  results  in  the  following  expression  relating  pressure 

(P)  to  height  (z): 

P  ,  Po  («.  -_  «>*? 
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where  y  is  the  dry  adiabatic  lapse  rate,  To  is  the  standard  atmospheric 
temperature,  w  is  the  specific  humidity  and  Po  is  the  standard  atmospheric 
pressure,  (w  =  .622  (e/P)  ==>  P  =  .622  e/w) 

An  expression  for  the  actual  vapor  pressure  (e)  is  obtained  by  using  the 
definition  of  specific  humidity  and  substituting  the  above  expression  in  for 
P.  The  lifting  condensation  level  is  obtained  when  the  following  two 
equations  are  equal  to  each  other: 

wPo   (To  -yz  )R 
e  =  '  '  ' 


.  622       T 


o 


es  =  6.11  exp  [  £-  (~-  - )] 

Rw  273   To  -yz 

After  a  series  of  substitutions,  using  Taylor  series  expansion  of 
In  x ~ x-1  and  several  algebraic  manipulations,  the  following  equation  is 
derived  giving  the  lifting  condensation  level: 

WPo  wPn 

zLCL   =  To    (102.041  -  14.6429   (23.3058  -  ln(—  ----)   -    [HK^gjjjg) 

oo  o™.       75418.2   ,l-  ^ 
-  23.30582 ]2    )) 

Only  surface     data  values     of  temperature     (To),    pressure     (Po )      and     specific 
humidity   (w)   are  required  to  find  the  LCL. 
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C.   User's  Guide 

1)   Turn  on  calculator  and  set  printer  to  the  NOR  mode. 


2) 


XEQ)  LCL 

3)  Prompt  reads,  SFC  TEMP  °C? 

4)  Type  in  a  surface  temperature. 


5)   Press  R/S  to  enter  the  temperature. 


6)  Prompt  is  SFC  PRES  MB? 

7)  Enter  a  surface  pressure. 


8)    R/S 


9)  Prompt  reads,  SP  HUMID  G/G? 

10)  Enter  a  specific  humidity  value. 


11)   R/S 


The  lifting  condensation  level  will  appear  with  its  value  in  meters.   If 
the  LCL  occurs  at  the  surface,  the  word  SURFACE  will  be  displayed. 

D.   Examples 

Examples  A,  B  and  C  in  Table  VI  illustrate  the  use  of  the  program  "LCL" 
in  estimating  the  lifting  condensation  level  from  a  surface  temperature, 
pressure  and  specific  humidity.  In  example  A,  a  surface  temperature  of  9°C, 
a  surface  pressure  of  1010  mb,  and  a  specific  humidity  of  .006  g/g  are  entered 
in  the  "LCL"  program.  The  data  result  in  a  lifting  condensation  level  of 
318.9m  which  is  verified  by  a  thermodynamic  diagram.  (See  Figure  10.)  In 
example  B,  the  lifting  condensation  level  occurs  at  the  surface  and  in  example 
C,  the  LCL  is  at  637.4  m. 
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Table  VI 
Examples  of  the  Lifting  Condensation  Level  Program. 


S^C  TEMP  C? 
3. 


.  0060 
LCI,  H,=3i3,8959 


SFC  TEMP  C? 

6.9988 
SFC  PRES  MB? 

1..  929. 9999 
SP  HUMID  G/G? 

.9849 
LCL,  H..  =637. 3732 
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B 


SFC  TEMP  C? 

.5800 


SrC  PRES  mo  SFC  PRES  MB? 

I;  819.9000  1..  998. 9900 

SP  HUMID  G/G?  SP  HijHir'  G/G' 


.9040 

SURFACE 


988m 


111m 


Figure  10 

Lifting  condensation  levels  determined  by  surface  temperature,   pressure 
and  specific  humidity. 
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E.   Program 


Table  VII 
Lifting  Condensation  Level  Program, 


Ml 

HBL  "LCL" 

92 

23.3953 

97 

?Tfl  25 

04 

*SFC  TE^P  C?" 

95 

PRO^P" 

96 

273. 15 

e? 

+ 

MX 

STO  21 

89 

"SFC  pR'ES  Ng-": 

19 

PROMPT 

1 1 

STO  22 

1? 

"SP  MiJH ! D  Q-G 

17 

PR0HFT 

i4 

STO  23 

15 

RCL  22 

\b 

* 

17 

3.39042 

iC 

i  'J 

19 

Ln 

■-•n 

STO  2* 

V' 

"■:■-•  1    -,C 

L  * 

F.VL   C  J 

22 

- 

"7 

>-*2 

75418.2 

V  ~ 

RCL  2: 

2  b 

00  c 


3  SURi 


29  s  T  ij  ££. 

30  CHS 

31  RCL  24 

32  CHS 

33  * 

3<*  RCL  25 

75  * 

36  14. 6*29 

37  * 

33  192.041 

39  XOY 

4f  - 

i)     Dpi     v' 

42  * 

43  STO  39 

44  0 

45  XOY 

46  X<=Yv 

47  ijJO   "SrP" 
43  "LCL,   H,=' 

49  PRCL  30 

50  hVIEh 

51  GTO  "END" 
52*LEL  "SFC" 


c- 


ippqc 


54  SVIEW 
55*L8L  "END" 
56  .END. 
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